Macroscopic description of particle systems with non-local density-dependent 
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In this paper we study macroscopic density equations in which the diffusion coefflcient depends 
on a weighted spatial average of the density itself. We show that large differences (not present 
in the local density-dependence case) appear between the density equations that are derived from 
different representations of the Langevin equation describing a system of interacting Brownian par- 
ticles. Linear stability analysis demonstrates that under some circumstances the density equation 
interpreted like Ito has pattern solutions, which never appear for the Ifanggi-Klimontovich inter- 
pretation, which is the other one typically appearing in the context of nonlinear diffusion processes. 
We also introduce a discrete-time microscopic model of particles that confirms the results obtained 
at the macroscopic density level. 
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In the mathematical description of many ecological 
systems the mobility of the species depends on the lo- 
cal density or concentration of the individuals [10. A 
standard example of this is given by models for animal 
or insect dispersal where it is usually considered that 
there is an increase of the diffusion coefficient due to 
population pressure. These models are written down as 
an extension of the diffusion equation which includes a 
diffusivity depending on the population density. Under 
the generic name of nonlinear diffusion equations similar 
equations arise in other contexts like in the flow through 
a porous medium, bacterial dynamics, and transport in 
plasmas 2, 3]. 

Generally, the nonlinear diffusion equations for the 
density correspond to a coarse-grained description of the 
more fundamental particle stochastic dynamics. Simi- 
larly to the case of diffusion in nonhomogeneous media 
where the diffusion coefficient varies in space, a heuris- 
tic derivation of the nonlinear diffusion density equations 
has to deal with the Ito-Stratonovich dilemma (To 
be precise one should speak about the interpretation of 
the Langevin equation since, as we will see later, the 
Stratonovich realization does not appear in our study and 
one has instead the Hanggi-Klimontovich interpretation 
of the Langevin equation 6].) In the context of noise- 
induced phase transitions this is of a major importance 
0, and this work pretends to show that it is also spe- 
cially relevant when one considers nonlocal interactions 
between the particles. 

In a biological setting, nonlocal effects have been 
widely considered since they account for the interaction 
between individuals that are separated in space |^] . More 
recently, there has appeared interest in models where the 
mobility depends on an averaged density rather than on 
its value at a point . The spatial non- locality is there 
introduced in order to take into account the finite-size 
of the particles or the senses (visual, hearing, etc..) of 



the organisms, that somehow determine their mobility. 
The idea behind is that the diffusivity of any particle 
may depend on the configuration of other particles in its 
vicinity. 

Thus, in this work we study a nonlocal version of a gen- 
eral density-dependent diffusion equation for the num- 
ber density of individuals. We focus on its pattern and 
clustering formation properties. Making linear stability 
analysis we first show that, contrarily to what happens 
in the local case, completely different results can be ob- 
tained for the Ito and Hanggi-Klimontovich descriptions. 
In fact, spatial patterns can emerge for the first but never 
for the second. Then, we introduce two microscopic in- 
teracting particles systems whose continuum evolutions 
are given by the Ito and Hanggi-Klimontovich density 
equations. Numerical simulations of these systems con- 
firm the calculations made with the continuum density 
equations. 

Let us consider a population dynamics model of rele- 
vance for insect dispersal. The diffusion coefficient de- 
pends on the local density of particles, i.e., D = D{p), 
where p(x, t) is the population density (we assume that 
D is not an explicit function of the spatiotemporal co- 
ordinates). We only study positive density-dependence, 
which is the most widespread hypothesis and indicates 
that competition increases dispersing because individ- 
uals have better survival rates by leaving high-density 
areas [13 ■ Mathematically, D'{p) = dD/dp > 0. 
The mean-field density descriptio_iis_iii _Jto and Hanggi- 
Klimontovich interpretations are 
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The first one corresponds to Ito and the second to 
Hanggi-Klimontovich, and in the following we will de- 
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note them, respectively, as IE and HKE. The HKE is 
the one that is usually considered when modeling at the 
macroscopic level P, , since one writes down the den- 
sity evolution as a flux equation dtp ~ — V • J with the 
density-dependent flux J = —D{p)Vp. But this is not a 
definitive argument and, as it is discussed in the IE 
should be the correct one for a continuum description of a 
discrete population dynamics model of non-overlapping 
generations. A proper derivation of density equations 
from the fundamental stochastic particle dynamics, once 
a right interpretation for this is assumed, can be con- 
sulted in 12]. Note that they are stochastic, but in our 
study we neglect the fluctuations by considering a kind of 
mean- field approximation [l^ . Let us assume an initial 
constant density p(x, 0) = po. We now perform a linear 
stability analysis of the stationary homogeneous (identi- 
cal to Po because particle number is conserved) solution 
for both IE and HKE. Writing p{x,t) = po + eilj{x,t), 
where e is a small parameter, and tp space-time de- 
pendent perturbation, we obtain (we label the equations 
with IE or HKE): 

(IE) dtij{-K,t) - (i5(po)+Po^'(po)) (3) 
{HKE) dt^j{-K,t) = i?(po)VV- (4) 

Note that in both the perturbation follows a simple dif- 
fusion equation and there is, therefore, no instability in 
the density equation that could give rise to the forma- 
tion of patterns or aggregates of particles. From this 
point of view, the Ito and Hanggi-Klimontovich frame- 
works (Eqs. and |2Il) arc equivalent. 

Now let us consider that the diffusivity depends on 
an averaged density over the whole system, yo(x, t) — 
J^a drG{x — r)p(r, t), where the kernel, G(x — r), repre- 
sents the effect of population density at r on the density 
at X. It embodies the forces of attraction and repulsion 
of neighbors 0|. We assume that G is normalized so 
that /sRd rfrG(r) ~ 1. Then let us consider D = D{p) 
in Eqs. (|1I2|I and make again a linear stability analysis 
around po for the now nonlocal HKE and IE, 

{HKE) 

efti^(x, t) = V [D{po + e J drG(x - r)V;(r, t)) 

V(po + e^)], (5) 

so that to first order in e the perturbation simply diffuses, 
dt^p — D{po)'V^ip, and no instability happens in HKE. 
However, when we consider the Ito description we have: 

{IE) e9tV(x, t) 

= [D {po + ej drG(x - r)V(r, t)) {po + eV(x, t))] 

= eV^ [D{po)ij + poD{po) J drG(x - r)^(r, t)] . (6) 

Considering a harmonic perturbation t) — exp(Af: + 
ik • x) it is not difficult to obtain the following dispersion 
relation 

A(fc) = -k^ (z?(po) + PoD'{po)G{k)) , (7) 



where k = |k| and G(k) = J drG{r)e~^^''^ is the Fourier 
transform of the kernel. The important fact with Eq. Q 
is that depending on the kernel function (note that the 
form of the diffusivity functional is almost irrelevant once 
we have assumed that it is an increasing function of the 
density, D'{p) > 0), the perturbation growth rate can 
be positive, giving rise to the aggregation of particles or 
formation of spatial patterns |3, llJI- As we have just 
seen, this is absolutely different to the result obtained in 
the Hanggi-Klimontovich framework, Eq. @ and below. 

Next we introduce a model, choosing the density and 
the kernel functions, with a clearly oriented ecological 
application and where A can take positive values. We 
study the distribution of organisms over a spatial area 
and thus we restrict to two spatial dimensions. A typical 
form for D{p) in animal dispersal models is Do{p/pq)p 
where Dq and p are positive real numbers 0,0 j and po is 
a reference density which we take equals to the initial one 
without loss of generality. Because of the normalization 
of the kernel function the averaged density, p, has the 
same dimensions as p. Concerning the kernel, we take 
the usual top-hat function 14]: G(r) = 1/{t:R'^) if |r| < 
R and G(r) = otherwise, which introduces a typical 
interaction length, R, in the system. One obtains G{k) = 
2Ji{kR)/{kR) so that the dispersion relation takes the 
form 

Ji is the first order Bessel function, and it is clear from 
Eq. (|SJ) that the value of Dq only sets the time scale of the 
system. The onset of pattern formation is A > which 
numerically is approached for p larger than pc ~ 7.6. In 
Fig. n we show a long-time density solution for the IE 
(left) and HKE (right) with this kernel and D functions, 
and p — 9. Note the hexagonal pattern formed in IE and 
the homogenous solution for HKE, i.e., the numerical 
results confirm the different spatial structures obtained 
for the two descriptions. We have checked that this is so 
for all p, as can be seen below in Fig. 13 

Up to now we have shown that, for nonlinear nonlo- 
cal diffusion phenomena, the different interpretations, in 
terms of density equations, of the same system dynamics 
lead to very distinct results. One may ask if this also 
happens at the level of the discrete interacting particle 
dynamics. Since the IE and HKE are the mean-field den- 
sity descriptions derived from the Langevin equation of a 
system of Brownian interacting particles, we expect that 
the same differences encountered before between both 
descriptions appear. We will show that this is in fact 
the case for two time-discrete microscopic systems which, 
in the continuous time limit, are the Ito and Hanggi- 
Klimontovich Langevin equations corresponding to the 
IE and HKE with the top-hat kernel and the typical dif- 
fusivity functional for insect dispersal. 

Therefore, the problem is to find a discrete dynani- 
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FIG. 1: Steady number density of particles with nonlo- 
cal density-dependence diffusivity. Left is Ito and right 
Hanggi-Klimontovich. The diffusivity functional is D{p) — 
Do{p/ poY , and we use the top-hat kernel with R = 0.1. The 
values of the other parameters, in both plots, are Dq = 0.0001, 
p = 9, po = 2000, and system sizes in x and y are equal to 1. 



FIG. 2: Spatial distribution of particles at large time. Left, 
is for Ito and right for Hanggi-Klimontovich. We use the 
same parameters as in Fig. i.e., Do = 0.0001, p ~ 9, and 
R — 0.1. The number of particles is = 2000, for both plots, 
and dt = 0.01. For the right panel we take a square grid (to 
calculate the gradients of the g function) with lattice spacing 
dx = 0.005. 



ics from which deriving HKE and IE. Following the dis- 
cussion m [13 (Section 5.4.2) we first take a population 
dynamics system of non-overlapping generations whose 
Langevin equation, in the time continuum limit, should 
be interpreted in Ito. The system consists in N parti- 
cles with positions ri{t) = {xi{t),yi{t)), i = 1,...,N, in 
a two-dimensional periodic box of size L — 1 evolving as 
follows 



r,{t + At) - r,{t) + ^/2Do{NR{r^{t))/N)PAt'n'^'\t), (9) 

where Dq and p are real positive numbers. At is the 
time step, Nii{ri{t)) is the number of particles at a dis- 
tance less than R of particle i at time t, and T7*^*'(i) = 
,riy^) is a Gaussian White noise with correlations 
(va\t)r]l''\s)'^ = SijSabSts- As already mentioned, in 
the limit At — > one obtains an Ito Langevin equation. 



-^ = V2Do{Nn{rm}/Nyv^'Ht), (10) 

and one can calculate 0| that the corresponding 
density equation of this Langevin is the nonlocal IE with 
the top-hat kernel and the density functional of insect 
dispersal. For the Hanggi-Klimontovich case, the first 
step is to interpret eq. (|l()|l in this sense so that it is 
equivalent to the following Ito Langevin equation jl5l | : 

dr. 



^=5(r.W)Vr..g(r,(t)) + v/2i^o(A^fl(r.(<))Wr,«(t), 

(11) 



where g{r,{t)) = ^2DoiNR{r,{t))/N)P. 

Thus, our discrete particle systems equivalent to the 
IE and HKE are given by Eq. ^ and the discrete-time 
version up to order At of Eq. Hll|l . respectively. Note 
that, since our simulations are off-lattice, i.e., particles 
move in the continuum space, we have to consider an 
auxilary grid to compute the gradients of g{ri{t)) for 



the Hanggi-Klimontovich case. In Fig. ||2Jl we show the 
long-time spatial distribution of particles for the same 
parameter values used in Fig. (Q. Left corresponds to 
Ito and right to Hanggi-Klimontovich. The particles fol- 
lowing the time-discrete-dynamics which give rise to the 
Ito Langevin equation arrange in an hexagonal pattern 
of clusters, while in the Hanggi-Klimontovich case the 
particles are randomly distributed in the space. This is 
completely in accord with what it happens in the macro- 
scopic description. 

A more quantitative comparison between the macro- 
scopic and microscopic dynamics follows. The character- 
ization of the pattern can be performed via the structure 
function, whose maxima identify periodic structures. In 
the case of the continuum density, the structure function 
in the steady state, S'c, is computed as the modulus of 
the Fourier transform of /9(x, t) averaged spherically and 
in time. In the case of the particle system the structure 
function, which is related but not identical to Sc, is cal- 
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culated as Sd{k) 



,{t) 



/N ), where Xj(t) is 



the position vector of particle j, k is a two-dimensional 
wave vector with modulus fc, and the average indicates a 
spherical average over the wave vectors with modulus fc, 
a temporal average in the steady state, and an average 
over many different realizations. In Fig. O we plot the 
maximum (for fc > 0) value of Sc (left panel) and of Sd 
(right) versus the parameter p. With circles we plot the 
computed ones for IE and with squares for HKE. It is 
clear from the plot that the value of p at which the tran- 
sition to pattern, in Ito, occurs is the same for both, and 
coincident with the numerical value obtained from stabil- 
ity analysis. Moreover, for any value of p no peak (i.e., 
no pattern) is ever observed in the Hanggi-Klimontovich 
prescription. 

In summary, density evolution equations with density- 
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FIG. 3: Maxima of the structure functions for different values 
of tlie parameter p. In the left panel we plot the maximum 
value of Sc, i.e., the corresponding to the continuum density 
equation. It is normalized with the number of points of the 
grid used to make the computation, 10201. Circles are for 
the Ito prescription and squares for the Hanggi-Klimontovich 
one. On the right, we plot, also vs p, the maximum value of 
Sd, i.e, of the structure function for a distribution of particles. 
Also circles are for Ito and squares for Hanggi-Klimontovich. 
The rest of the parameters for both plots are Do = 0.0001 
and R = 0.1. For the Hanggi-Klimontovich particle dynamics 
we use a lattice with dx = 0.005. 



dependent diffusion coefRcient are standards in the 
macroscopic modeling of biological, physical and chem- 
ical processes. Typically, they are postulated as 
9tp(x,t) = —V • J with J the population flux given by 
3 = —D{p)'V p. In this way one obtains a macroscopic 
density equation in the so-called Hanggi-Klimontovich 
framework, which is different to the Ito one, Eq. 
Though this last expression rarely appears in the liter- 
ature, this is no reason to exclude it 0, and one has 
to resort to the real microscopic dynamics to infer the 
proper density equation describing the system dynam- 
ics. In fact, a proper macroscopic equation for popu- 
lation dynamics models of non-overlapping generations 
goes through Ito this work we have shown that, in 

the context of biological population dynamics but easily 
extended to others, no relevant differences, at the level 
of clustering or pattern formation, arise between both 
equations when the functional D{p) is local. However, 
when D is nonlocal the differences are rather important, 
and non-homogenous spatial structures can appear in 
Ito but never in Hanggi-Klimontovich. Furthermore, we 
have discussed the discrete dynamics of non-locally inter- 
acting particle systems from which the Ito and Hanggi- 
Klimontovich density equations are derived, confirming 
also at this level the large differences encountered in the 
macroscopic description. It is thus clear that the main 
message of this work is that it is specially relevant in the 
case of systems with nonlocal density-dependent diffusiv- 
ity that the macroscopic description should be properly 
derived from the discrete particle model. The inclusion 
of reaction terms and further extensions of the model will 



be pursued in the future. 
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